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£ — , | Low energy super symmetric models provide a solution to the hierarchy problem 

and also have the necessary ingredients to solve two of the most outstanding issues in 
cosmology: the origin of dark matter and baryonic matter. One of the most attractive 
features of this framework is that the relevant physical processes are related to inter- 
actions at the weak scale and therefore may be tested in collider experiments in the 



near future. This is true for the Minimal Super symmetric Standard Model (MSSM) as 
well as for its extension with the addition of one singlet chiral superfield, the so-called 
nMSSM. It has been recently shown that within the nMSSM an elegant solution to 
both the problem of baryogenesis and dark matter may be found, that relies mostly 
on the mixing of the singlet sector with the Higgs sector of the theory. In this work 
we review the nMSSM model constraints from cosmology and present the associated 
collider phenomenology at the LHC and the ILC. We show that the ILC will efficiently 
probe the neutralino, chargino and Higgs sectors, allowing to confront cosmological 
observations with computations based on collider measurements. We also investigate 
the prospects for a direct detection of dark matter and the constraints imposed by the 
current bounds of the electron electric dipole moment in this model. 



1 Introduction 



In spite of the excellent agreement of the predictions of the Standard Model (SM) with ex- 
perimental observables, there is still a strong physical motivation for the presence of new 
physics at the weak scale. The main reason is the belief that the Higgs mechanism of the SM 
is only an effective description of a more fundamental mechanism of electroweak symmetry 
breaking, in which the origin and stability of the electroweak scale must be explained. Su- 
persymmetric extensions of the SM [1,2] allow such a mechanism. In the simplest extensions 
of the SM, the theory remains perturbative up to scales of the order of the Planck scale 
and the weak scale is stable under quantum corrections at all orders in perturbation theory. 
Moreover, the electroweak symmetry is broken radiatively, providing a correlation between 
the weak scale and the soft supersymmetry breaking scale. 

Another strong motivation is a solution to the problem of dark matter. If a discrete 
symmetry, R-Parity, is imposed in the theory, the lightest supersymmetric particle (LSP) is 
stable. In the simplest models, it is also neutral and weakly interacting, with an annihilation 
cross section of the order of the one necessary for the LSP to become a good dark matter 
candidate. 

It has been also realized that supersymmetry may also lead to a solution of another 
outstanding cosmological problem, namely, the baryon-antibaryon asymmetry [3,4]. In the 
MSSM, a solution to this problem by weak scale physics demands a light stop as well as 
non-vanishing CP-violating phases in the chargino-neutralino sector. The presence of a 
light stop, with mass smaller than the top quark mass, tends to push the lightest CP-even 
Higgs mass to values close to the ones restricted by experimental searches. Therefore, this 
scenario is highly constrained and it will be probed at the Tevatron and the LHC colliders 
in the near future, as well as in electric dipole moment and direct dark matter searches [5,6] 
and the ILC [7]. Supersymmetry also plays a relevant role for the generation of the baryon 
asymmetry in the so-called soft leptogenesis scenario, which again demands very specific 
values for the soft supersymmetry breaking parameters in the neutrino sector [8], which can 
be naturally obtained only in certain supersymmetry breaking scenarios [9]. 

A less constrained solution to the baryon asymmetry may be found within extensions of 
the MSSM with additional chiral singlet superfields. In particular, it has been recently shown 
that a particularly simple extension of the MSSM, the so-called nMSSM [10-14] includes all 
the relevant properties to lead to a successful generation of the baryon asymmetry at the 
weak scale [15,16], without the need of light squarks and with values of the lightest CP-even 
Higgs mass which are naturally above the current experimental bound [17]. 

One interesting feature of the nMSSM is that the lightest neutralino is usually an ad- 
mixture of the fermion component of the singlet field and the Higgsinos. Due to the absence 
of explicit mass terms in the Higgsino-singlino sector, its mass is bounded to be below 70 
GeV. A relaxation of this upper bound may only be obtained by allowing the trilinear Higgs- 
singlet coupling to become strong at scales lower than the Grand Unification (GUT) scale, 
as happens in the so-called fat Higgs models and gauge extensions* of the nMSSM [18]. In 

*Such gauge extensions open new possibilities for baryogenesis [19], but this possibility will not be pursued 
any further in this paper. 



1 



this work, we shall demand perturbative consistency of the theory up to the GUT scale, and 
therefore the aforementioned bound applies. For light neutralinos, the annihilation cross 
section is dominated by s-channel Z° boson exchange and a neutralino relic density consis- 
tent with all experimental bounds, may only be obtained for values of the neutralino mass 
of about 30-45 GeV. 

Such a light neutralino has relevant phenomenological consequences. For instance, the 
lightest CP-even Higgs tends to decay predominantly into neutralinos, leading to invisible 
signatures at hadron and lepton colliders [16]. Moreover, one of the predictions of this model 
is the presence of relatively light charginos and neutralinos with masses smaller than about 
200 GeV. Detection of these weakly interacting particles at hadron colliders becomes difficult 
if there are no strongly interacting superpartners with masses below the TeV scale. Assuming 
that the relevant CP-violating phases are in the gaugino masses, successful baryogenesis may 
be easily achieved for weak-scale gaugino masses of the order of 100-300 GeV, and at least 
one relatively light squark in order to avoid the so-called Giudice-Shaposhnikov suppression 
of chiral charges [20]. We will assume the presence of top and bottom squarks with masses 
of order of 500 GeV, and the presence of a gluino with mass dictated by the gaugino mass 
unification condition. The phenomenological properties at hadron colliders depend strongly 
on these last assumptions. 

In addition to discussing the prospects for discovery of the light Higgs and superpartner 
states at the LHC, we will perform a detailed phenomenological analysis of nMSSM physics 
at the ILC. In this case the values of the strongly interacting particle masses become less 
important. The same is true for the analysis of the direct dark matter detection properties 
and the constraints coming from the bounds on the electron electric dipole moments that 
we investigate in this work. 

The article is organized as follows. In section 2 we present an overview of the nMSSM and 
present experimental constraints to define a benchmark scenario for the model. In section 
3 we analyze the corresponding collider phenomenology at the LHC and the ILC, which 
we then connect with the cosmology constraints related to the dark matter relic density 
and baryogenesis in section 4. In particular, we present the constraints coming from direct 
searches of dark matter and the non-observation of an electron electric dipole moment. We 
reserve section 5 for our conclusions. 

2 General properties of and constraints on the nMSSM 
2.1 Overview 

One of the original motivations for a singlet extension of the MSSM is the difficulty to explain 
why the /i-parameter in the superpotential is of order of the electroweak scale, instead of the 
GUT or Planck scale. This problem is circumvented by introducing a singlet chiral superfield 
S, and generating a weak-scale /z-parameter through the vacuum expectation value of the 
scalar component S of the singlet. This requires the inclusion of a triple-Higgs coupling term 

W x = XSHi ■ H 2 (1) 
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in the superpotential, where H 12 are the superfields of the two Higgs doublets, 
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This triple-Higgs coupling also helps to push up the mass of the CP-even Higgs boson respon- 
sible for electroweak symmetry breaking, so that the bound from LEP, 
m h o > 114.4 GeV [17], can be satisfied without substantial fine tuning. Finally, the new 
contribution increases the strength of the electroweak phase transition, thus allowing the 
possibility of electroweak baryogenesis in a large region of parameter space [16]. 

In a general singlet extension, the superpotential can also contain a triple-self-coupling, 
a mass term and a tadpole term for S, 



where y u ,d,e are the 3x3 Yukawa coupling matrices. 

In the Next-to-Minimal Supersymmetric Standard Model (NMSSM), a discrete Z 3 sym- 
metry is imposed, which forbids the tadpole and mass terms and only leaves dimensionless 
couplings in the superpotential. However, once the singlet scalar acquires a vacuum ex- 
pectation value, the Z 3 symmetry is broken and unacceptable domain walls can be gener- 
ated [10,21]. 

In the Nearly Minimal Supersymmetric Standard Model (nMSSM) [11-14], a Z 5 or Z 7 
R-symmetry is introduced in the superpotential and the Kahler potential, which forbids 
quadratic and cubic self-interactions of the singlet. However, at higher loop orders, a tadpole 
term is induced due to supergravity effects. It is suppressed at the six (Z5) or seven (Z7) loop 
level, so that it is naturally of weak-scale order and does not destabilize the hierarchy [12,13]. 
The discrete symmetries also ensure that rapid proton decay does not occur, while Majorana 
neutrino masses are allowed. Furthermore the lightest supersymmetric particle, usually 
the lightest neutralino, is quasi-stable on cosmological time-scales [16], even without the 
introduction of R-Parity. The superpotential of the nMSSM is thus 



In the following, the phenomenological properties of the nMSSM will be analyzed in more 
detail. 



In the nMSSM superpotential, the parameters A, mi 2 and yt can in general be complex. 
However, the phase of A does not generate physical CP-violating effects, since it can be 



-s -s -s -s -s III /V 

WWm = XSHi ■ H 2 + kS 3 + m N S 2 + -^S 

A 

+ y u Q-H 2 U c + y d Q-H 1 D c + y c L-H 1 E c , 



(3) 



H^nMSSM = XSH 1 ■ H 2 H — 

+ yuQ-H 2 U c + y d Q-H 1 t) c + y c L-H 1 E c . 



(4) 



2.2 CP violation 
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Table 1 : Charges of superfields under the U(l ) R and U(l ) P q symmetries of the superpo- 
tential. Note that the U(l ) R charges of the fermionic components differ by 1 from those of 
the superfields. 



made real by suitable gauge rotations [22]. As in the Standard Model, the Yukawa couplings 
yi lead to one physical phase in the CKM quark mixing matrix, which is constrained by 
present data from many heavy-flavor experiments. The phase of m\ 2 will be addressed 
below. 

Beyond the superpotential, new complex phases can appear through supersymmetry 
breaking. The soft supersymmetry breaking Lagrangian reads 



£ soft = mjHlHi + m\H\H 2 + m 2 s \S\ 2 + (t s S + a A S#i • H 2 + h.c.) 
+ {M-JSB + M 2 W -W + M 3 gg + h.c.) 
+ m 2 ~ql- q L + m\\u K \ 2 + mj~|(i R | 2 +m~P L - l L + m?|e R | 2 
+ (VuA u q L ■ H 2 u* K + y d A d q L ■ H 1 d* R + y c A c l L ■ H 1 e* R + h.c). 



(5) 



Here Hi, S, qi,, Mr, ^r, II, eR are the scalar components of the superfields Hi, S, Q, U, D, L, E, 
where the quark and lepton fields exist in three generations (the generation index has been 
suppressed in the formula). B,W,g denote the fermionic components of the gauge super- 
multiplets. Among the soft breaking parameters, a\, t s , Mi i2j3 and A u ^ c can be complex. 
However not all their phases are physical. To see this, one can observe that the superpotential 
is invariant under an U(1)r symmetry, with the charges listed in Tab. 1. In addition, it obeys 
an approximate Peccei-Quinn symmetry U(1)pq, which is broken by the singlet tadpole term 
oc m\ 2 . Both U(1)r and U(1)pq are broken by some of the supersymmetry breaking terms. 

With the help of the U(1)r and U(1)pq, the fields can be rotated so that the phases of 
two parameters become zero. By analyzing the charges, it can be seen that the following 
products remain invariant under both R- and PQ-transformations: 

arg(m 2 2 t*a A ), 

arg(m 2 2 t s *M,), i = 1,2,3, 

arg(m 2 2 t*A u ), (3 generations), (6) 

arg(m 2 2 t*Aj), (3 generations), 

arg(m 2 2 t*A e ), (3 generations), 

corresponding to 13 physical CP- violating phases in addition to the CKM phase. Without 
loss of generality, the phases of m\ 2 and t s can be chosen real, so that the physical phases 
are transferred into a\, Mi j2 ,3 and A u ^ c . 
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In this work, for simplicity, gaugino unification is assumed, so that M 1 : M 2 : M 3 I : 
2:6. In this case, the gaugino masses carry one common phase, 0Mi = 0m 2 — 0m 3 = 0m- 
To simplify the analysis further, the phases in Ai,d,e and a\ are set to zero. 



2.3 Higgs sector 

The scalar potential at tree-level receives contributions from F-, D-terms and soft super- 
symmetry breaking. To avoid charge and color breaking vacuum solutions and unacceptably 
large neutrino masses, the squarks and sleptons must not receive vacuum expectation values. 
In this case, the Higgs potential at tree-level reads 

V = V F + V D + V soft , (7) 



Vv = 



XHi ■ H 2 + 



m 2 2 



A 



2 

+ \XS\ 2 (HlH 1 + H 2 l H 2 ), (8) 



Vd = ^(HlH, - H\H 2 f + g ^\H\H 2 \\ (9) 
V; oft = m\H\Hi + m\H\E 2 + m 2 s \S\ 2 + (t s S + axSHi ■ H 2 + h.c), (10) 

where g 2 = g 2 + g' 2 , and g, g' are the SU(2), U(l) gauge couplings, respectively. As men- 
tioned in the previous section, we constrain all parameters in the Higgs sector to be real. 
Nevertheless, we allow complex phases in the gaugino sector, which in turn would generate 
small CP-violating phases in the Higgs sector through radiative corrections. However, the ef- 
fect of gaugino loop contributions to the Higgs effective potential is generally sub-dominant, 
compared to the top-stop loops, and thus can be neglected in the following discussion. 

In the zero-temperature vacuum state, the neutral Higgs components acquire non-zero 
vacuum expectation values, 

(S)=v„ (H° 1 )=v 1 , (H° 2 )=v 2 . (11) 

It is useful to define tan/3 = v 2 jv\. The minimization conditions for electroweak symmetry 
breaking give at Born level 

— 2 

mj = -(m? 2 + a x v s )^- - 9 -{v\ - v\) - \ 2 {v\ + v 2 ), (12) 

— 2 

m\ = -{m\ 2 + a x v s )^ - 9 -{v\ - v\) - \\v\ + v 2 ), (13) 

2 V 1 V 2 t s 2 2 

m s = -a x A v . (14) 

The physical mass eigenstates of the Higgs sector consist of three neutral CP-even scalars 
5*1,2,3, two neutral CP-odd scalars P ij2 and one charged scalar H ± . Both the CP-even and 
CP-odd neutral scalars share a component of the singlet scalar S. 

In the sector of the CP-odd scalars, one combination of the two doublets form the Gold- 
stone mode G°, which is absorbed into the longitudinal mode of the Z boson. The other 
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linear combination of the two doublets, A , mixes with the imaginary part of the singlet S 
to give the two physical states. In the basis (A , 5smS), their mass matrix reads 



Mi 



Ml 
-a x v 




where 



Ml = (m\ 2 + a x v s ), 

$2/3 



(15) 



(16) 



and s 2 [3 = sin 2/3. The mass eigenstates Pi 2 are related to this basis by the mixing matrix 
O p , 

'Pi\ „J A 



P'2 I \ ^sm S 



(17) 



For large values of Ma, corresponding to large negative values of m\ 2 , the mixing between 
the two CP-odd scalars is relatively small, with the light pseudo-scalar Pi being almost 
completely singlet, while the heavy pseudo-scalar P 2 has a very small singlet component and 
a mass of about m P2 ~ Ma.. These properties have strong consequences for Higgs searches 
at LHC and ILC, as discussed in the next section. 

Similarly to the CP-odd sector, the CP-even Higgs eigenstates are related to the doublet 
and singlet components by the mixing matrix s , 



(18) 



For Ma <C Mz, there are one heavy mass eigenstate S3 with small singlet component, and 
two light eigenstates Si, £2 with sizable singlet and doublet admixtures. 

The Higgs potential receives large radiative corrections, with the dominant contribution 
stemming from top-stop loops. These effects have been calculated in Ref. [13,23] and are 
included in this work. 

The influence of the parameter A gives an additional positive contribution to the mass 
of the lightest Higgs boson, so that the bound of about 114 GeV from LEP can be easily 
avoided even for moderate values of tan (5 and stop mass eigenvalues m t ~ , m t ~ . From the 
vacuum expectation value of S, v s , an effective /z-parameter is generated, \i = —Xv s . 



(s 1 \ 




( H°\ 


S2 


= O s 


HI 


\SsJ 




\$te Sj 



2.4 Chargino and neutralino sector 

The presence of the fermion component of the singlet superfield, the singlino S, leads to a 
fifth neutralino state. The chargino mass matrix in the basis (W ± ,H ± ) is 

M# = ( ^ , (19) 

x \V2cpM w -\v s 
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while the neutralino mass matrix in the basis (B°, W°, H®, H®, S) is given by 



/ M 1 




M 2 

C/3C W M Z 
-Sf3C W M z 







-Cfss w M z 
cpc w M z 



s p s^M z \ 
-S/3C W M Z 

\v s \v 2 

Aui 

Aui y 



M^o = 









\v 2 



(20) 



V 



o 



where the abbreviations sp = sin/5, eg = cos/?, % = sin# w , c w = cos#w have been used. 
The neutralino mass matrix is diagonalized by a unitary matrix N, such that 



Perturbativity of the trilinear coupling A up to the GUT scale requires that A < 0.8 [16]. 
For values of A < 0.8, the lightest neutralino has a large singlino component, and a mass 
below about 60 GeV. 

2.5 Experimental and astrophysical constraints 

The chargino and neutralino spectrum is constrained by new physics searches at LEP. For 
the lightest chargino, the LEP measurements lead to the constraint m-± > 104 GeV or 
a[e + e~ — > xfxf] x BR < 10 fb, where a is the production cross section at ^fs = 208 GeV, 
and BR the branching ratio into hadron or leptons. Similarly, the corresponding requirement 
for neutralinos is (m^o + m^o) > 208 GeV or <x[e + e~ — > X1X2] x BR < 10 fb [24]. Moreover, 
if the lightest neutralino is so light that m^o < M z /2, the measurement of the total Z-boson 
width imposes the limit BR[Z -> XiXi] < 0-8 x 10~ 3 [25]. On the other hand, if the nMSSM 
is to account for the dark matter relic density in the universe, the lightest neutralino mass 
is bound to be in the range 25 GeV < m^o < 40 GeV [16]. 

For the mechanism of electroweak baryogenesis to work successfully, the electroweak 
phase transition needs to be strongly first order, and new sources of CP-violation beyond 
the CKM matrix must be present. The electroweak phase transition can be made first order 
through a relatively large cubic term in the effective Higgs potential. In the nMSSM, similar 
to other singlet extensions of the MSSM, this cubic term is realized at tree-level by the a\ soft 
breaking parameter. As a consequence, a strongly first order phase transition can be achieved 
without relying on radiative contributions from a light stop as in the MSSM [16,26,27]. As 
a result, neither of the stop masses is constrained to be small ( < 150 GeV in the MSSM), 
but rather can amount to several hundred GeV. 

CP-violation in baryon-number violating processes can be introduced via complex cou- 
plings in the Higgs sector or through complex parameters in currents that couple sufficiently 
strongly to the Higgs potential. In the nMSSM, four physical phases can thus contribute to 
the baryon asymmetry, the phases of M ± , M 2 , A t and t s [27]. In this work, as emphasized 
before, it is assumed that the only non-zero phase (besides the CKM matrix) is a common 
phase of the gaugino masses, 4>m 1 = <Pm 2 — 0m 3 = 0m- 



N*M^oN' = diag(m^o, m^o, m^o, m^o, m^o). 



(21) 
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In the presence of CP-violating phases, the nMSSM parameter space is constrained by 
the experimental limits on the electric dipole moment (EDM) of the electron, neutron and 
199 Hg nucleus [28,29]. The upper bound on the electron EDM is derived from limits of the 
EDM of the 205 T1 atom. For the phases considered in this work, without mixing between 
the CP-even and CP-odd Higgs states, the CP-odd electron-neutron operator studied in [30] 
vanishes, and the 205 T1 EDM is due almost entirely to the electron EDM. This translates 
into a limit on the electron EDM of [28] 

| del < 1-9 x KT 27 e cm, (22) 

at 95% CL. 

The electron EDM receives potentially large contributions from one-loop sfermion-gaugino 
diagrams [31], which become small for large masses of several TeV for the first two genera- 
tion sleptons. A similar role is played by the squarks in the EDMs of the neutron and Hg 
atom. Two-loop contributions from diagrams with charginos and charged Higgs propagators 
can also be sizable, but just as in the MSSM [32] they can be suppressed below the current 
limits for sizable values of the pseudo-scalar masses and tan (3 ~ 0(1). However, since the 
two-loop contributions are enhanced by the same complex phase that generates the baryon 
asymmetry, the electron EDM bound presents a particularly severe constraint on electroweak 
baryogenesis. 

The one- and two-loop contributions to the electron EDM were calculated following 
Refs. [32] and [30,33], respectively. In the calculation of the one- and two-loop electron 
EDM contributions, we have included the Higgs and slepton mixing. 

Note that while slepton mixing is suppressed with the very small electron mass, the 
overall contribution to the electron EDM requires a spin flip and is thus also proportional 
the electron mass. As a result, slepton mixing can lead to non-negligible effects when the 
slepton masses are moderate. On the other hand, the contribution from the heaviest of the 
Higgs states S3 and P 2 is suppressed due to their large masses. 

For first and second generation sfermion masses of about 10 TeV, pseudo-scalar mass Ma 
higher than 500 GeV, and values of tan/5 ~ 0(1), which is in accordance with baryogenesis 
constraints, we find that phases 0m ~ are allowed by the electron and neutron EDM 
constraints. Since a sufficiently large baryon asymmetry could be generated for phases in 
the range 0m ~ 0.1 . . . 0.3 [16]*, this leaves a large window for this scenario to be realized in 
the nMSSM. 

Here we do not attempt to provide a theoretical model that can explain this pattern for 
the superpartner masses and parameters, but it has been pointed out in the literature that 
a roughly similar pattern might emerge within the framework of split supersymmetry [38]. 

t Different methods have been used to compute the baryon asymmetry in the MSSM [34-37]. The values 
derived in Ref. [16] are based on the method of Ref. [34]. While other methods [36] tend to lead to a 
higher baryon asymmetry, a more recent calculation leads to a lower one [37]. Furthermore, there are 
several approximations performed in the computation of the baryon asymmetry, leading to an uncertainty, 
of order one, on the necessary CP-violation phases required for the realization of this scenario. In section 
4.3.2, we provide a more general discussion of the possible sources of CP-violation, and how they affect the 
phenomenology of this model. 
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tan p 


A 






,1/3 


M A 


\Mi\ 


|M 2 | 


|M 3 | 


0M 


1.7 


0.619 


-384 


373 


157 


923 


122.5 


245 


730 


0.14 




m[/i,2 




m Llt2 


m E i,2 


™Q3 


m U3 


m D3 




m E 3 A,b,r 


10000 


10000 


10000 


10000 


10000 


500 


500 


500 


500 


500 -100 



Table 2: Parameters for the reference point A [16] used in this study. All dimensionful 
parameters are given in GeV. 

2.6 Benchmark scenario A 

In the following sections, the collider phenomenology and dark matter detection expectation 
will be analyzed in detail for the concrete parameter point A from Ref. [16]. The parameters 
of this scenario are summarized in Tab. 2. To avoid EDM constraints, the sleptons and 
squarks of the first two generations have masses of 10 TeV in benchmark scenario A. In 
contrast, the masses of the third generation sfermions are assumed to be around 500 GeV. 
There are no strong bounds for these masses, but stop masses with a few hundred GeV 
are favored by Higgs mass naturalness and baryogenesis. The masses and decay modes of 
the neutralinos, charginos, Higgs scalars, third generation squarks and the gluino for this 
scenario are listed in Tab. 3, Tab. 4 and Tab. 5. 

The relatively light neutralino and chargino spectrum is typical for a nMSSM scenario 
that is in agreement with electroweak baryogenesis and the relic dark matter density [16]. 
The lightest neutralino has a large singlino component, so that the bound from the invisible 
Z-boson width is evaded although m^o < M%/2. Nevertheless, the dark matter annihilation 
cross-section is dominated by s-channel Z exchange, since the combined mass of two Xi is 
close to the Z resonance, 2m^o ~ Mz- As the sfermions are relatively heavy, the other 
neutralino and chargino states mainly decay through gauge bosons. 

Owing the small mass of the lightest neutralino Xi, the light Higgs scalars Si, S 2 and Pi 
decay predominantly into two Xi- F° r the light CP-odd scalar Pi, which is an almost pure 
singlino state, this is essentially the only allowed channel, with all other contributions being 
far below 1% and thus negligible. 

Since the most relevant contribution to the annihilation cross section is coming from the Z 
pole, in Ref. [16] a simple calculation sufficed to determine the relic abundance of neutralinos 
for the benchmark cases.-'- Due to the precision nature that the ILC lends to our present 
work, we paid careful attention to reproducing, and improving, the earlier results. In our 
work, we include an improved superpartner spectrum calculation, all possible co-annihilation 
channels, all SM final states, and solve the relevant Boltzmann equation numerically for the 
relic density calculation. 

While our numbers closely agree with those of Ref. [16] there are some small deviations. 
Notably, our relic density calculation results in slightly higher values. When using the input 
values for benchmark A, we obtain a neutralino relic density of flh 2 = 0.131. The deviation 
of this from the latest WMAP central value Vth 2 = 0.111 [39] is comparable to the other 

■'•We thank A. Menon and D. Morrissey to provide us with their original codes. 
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Sparticle 



X? 

xl 
xl 



xl 



xl 



xt 
xt 



Mass m [GeV] 



33.3 

106.6 

181.5 



278.0 



324.4 



165.0 
319.5 



Width T [GeV] 



0.00004 
0.09 



1.5 



2.1 



0.136 
2.0 



Decay modes 



Xl 

xl 



xl 



xl 



zx\ 

SiX°i 

Pixl 

zxi 

zxl 
zxl 

W ± xl 

Six , 

Sixl 

S2XI 

PiX°i 

Pixi 
zxi 

zxl 
zxl 

W ± xf 

SiX? 
Sixl 
Sixl 

S2XI 
S2XI 

Pixl 
Pixl 
Pixl 



XI 

xt 



w + xi 
w+xi 

W + xl 

w + xl 

zxi 

Si Xi 
PiXt 



100% 
74% 
26% 
0.4% 

11% 

22% 

1% 

43% 

7% 

0.2% 

8% 

7% 

0.7% 

30% 

1.5% 

0.15% 

57% 

0.01% 

0.02% 

5% 

1% 

4% 

0.4% 

0.7% 

0.06% 



100% 

32% 

1% 

34% 

29% 

5% 

0.3% 



Table 3: Masses, widths and main branching ratios of the neutralino and chargino states 
at Born level for the reference point A (Tab. 2). 
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Sparticle 



Si 

s 2 



Pi 



Mass m [GeV] 



115.2 



156.6 



133.7 



Width T [GeV] 



0.044 



0.060 



0.008 



Decay modes 



Si 
S 2 



bb 

bb 

w + w~ 
zz 

xix\ 



Pi 



xixi 



8% 
92% 

2% 
17% 

2.5% 
69% 
10% 



100% 



Table 4: Masses, widths and main branching ratios of the light Higgs states, including 
one-loop top-stop corrections, for the reference point A (Tab. 2). 



theoretical and experimental uncertainties (higher order corrections to sparticle masses and 
annihilation cross sections, systematic errors, etc.) entering into this calculation. 

To demonstrate the uncertainty in the relic density calculation, we note that a variation 
of the neutralino mass by less than 1 GeV is enough to change the relic density by 0.025 for 
scenario A. Since the main focus of our work is collider phenomenology and a 1 GeV shift 
in the neutralino mass affects collider phenomenology negligibly, our results remain valid for 
all values close to VLh 2 = 0.111. 

Therefore, in spite of the discrepancy in flh 2 , we will stick to the scenario A in our 
numerical analysis. 



3 Collider measurements in the nMSSM 

The nMSSM can play a crucial role in baryogenesis and dark matter generation, depending 
on the parameters of the Higgs, chargino and neutralino sectors. In this section it is studied 
how the relevant particles can be discovered and the relevant parameters be determined at 
future colliders. 



3.1 The nMSSM at LHC 
3.1.1 Higgs physics 

In the benchmark scenario A, the light Higgs states Si, S 2 and P 1 mainly decay invisibly 
into the lightest neutralinos. This behavior is rather typical for many singlet-extended 
supersymmetry models [40]. A CP-even invisible Higgs boson can be discovered at the LHC 
through H^-boson fusion [41] and through associated Zh production [42]. The cross-sections 
and final state distributions for these processes in the nMSSM are the same as for the 
Standard Model Higgs boson, except for modified Higgs- Z couplings, that are given by 

Gzzs^G^ispOi + cpOi). (23) 
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Sparticle 



ti 



bi 



Mass m [GeV] 



522.1 



535.3 



498.8 



503.3 



730 



Width T [GeV] 



43.1 



29.8 



6.6 



12.5 



6.6 



Decay modes 



m 

txl 
txl 
txl 
txl 



bx 2 

txl 
txl 
txl 



tx 
txl 

bxt 
bxt 



bxl 
bxl 
bxl 
bxl 



tXi 

tX.2 



9 



12% 

7% 

11% 

20% 

1% 



bxt 47% 
: + 2% 
14% 
2% 
11% 
9% 
20% 
6% 
37% 



6i -> 6x? 0.5% 



12% 
12% 
0.1% 

5% 



-**r 57% 

- *X 2 ~ 15% 

6 2 -> 6x? 0.6% 

-6xS 4% 

-&X§ 9% 

-+6x2 0.2% 
7% 



56% 
23% 



661 35% 

66 2 34% 
tti 18% 
tt 2 13% 



Table 5: Masses, widths and main branching ratios of the third generation squarks and the 
gluino at Born level for the reference point A (Tab. 2). 
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pp — > 


99 hb\ 


b 1 b* 2 /b 2 b* 1 


6 2 &2 






t 2 t* 2 






<7[fb] 


2162 444 


3 


421 


357 


1 


312 


141 


138 



Table 6: Tree-level production cross-sections for strongly interacting super symmetric parti- 
cles at the LHC for the reference point A, see Tab. 2. The QCD scale is taken at the average 
mass of the produced particles, Q = (mj + m,y)/2. 

For the reference point A (Tab. 2), the couplings for the three CP-even Higgs states amount 
to 

Gzzsi — 0-87 x ^ZZH> 

G Z zs 2 = 0.49 x G™ u , (24) 
Gzzs 3 = 0.002 x G^h- 

Since both the Si and S 2 states have 0(1) couplings to the Z boson and are relatively light, 
with masses below 200 GeV, they can be produced with sizable rates. Based on the analysis 
of Ref. [41], it can be estimated that with only a few fb -1 of integrated luminosity, a 5er 
discovery in the IT-boson fusion channel can be achieved in the scenario A. 

For the case of a single Standard-Model-like Higgs boson, the mass of the Higgs can be 
determined from the ratio of the H^-boson fusion and Zh production rates [43]. On the other 
hand, in scenario A, the invisible Higgs signal receives contributions from both the S\ and 
S 2 . However, the observables in both IT^-boson fusion and Zh production are not sensitive 
to discriminate between one and two invisible Higgs states. As a result, from an invisible 
Higgs signal at the LHC, it is not possible to obtain information about the number and the 
masses of the CP-even Higgs bosons in the nMSSM. 



3.1.2 Supersymmetric particles 

Since the partners of quarks and the gluon, squarks and gluino, couple with the strong QCD 
coupling, they are produced with large cross-sections at the LHC. Charginos and neutralino 
can be generated in the decay cascades of squarks and the gluino with sizeable rates. In 
principle, charginos and neutralino are also produced directly in electroweak processes, but 
the cross-sections for these channels are small. Therefore in the following only squarks and 
gluinos are considered as primary supersymmetric particles. At the Tevatron, the typically 
large slepton masses necessary to suppress the electron electric dipole moment in this model 
lead to a reduced branching ratio for the decay of neutralinos and charginos into lepton final 
states. This makes their searches in the tri-lepton channel quite difficult, particularly for 
masses of the chargino and second lightest neutralino larger than 150 GeV, as are typical in 
the nMSSM scenario under analysis [44]. 

The relevant leading-order production cross-sections for squarks and gluinos at the LHC 
are summarized in Tab. 6, calculated with CompHEP 4.4 [45]. The analysis of x 2 produc- 
tion is experimentally particularly promising [46]. The neutralino x 2 * s produced in various 
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squark and gluino decay cascades, leading to a total cross-section for production with 
leptonic \\ decays of 30 fb. Here the most important channel is 

PP ^99, 9 -> bb* or bb -> bbx° 2 , (25) 

but direct production of sbottoms and stops via pp — > bb*,tt* also plays a role. According to 
Ref. [46,47], background from SM gauge bosons can be reduced by cuts on missing transverse 
energy and missing mass: 

• At least three jets with transverse momentum pf * > 150, 100, 50 GeV. 

• Missing energy # > max(100 GeV,0.2M cfr ) with M cS = $ + ELi^J- 

• Two isolated leptons with p\ cp > 20, 10 GeV. 

The remaining tt background is removed by subtracting events with two different-flavor 
leptons from events with same-flavor leptons. This procedure makes use of the fact that the 
tt background produces the same number of same-flavor and different-flavor lepton pairs, 
while the neutralino signal has only same-flavor lepton pairs. After these cuts practically 
no SM background is left, while a signal efficiency for x 2 production of better than 20% is 
achieved [46,47]. This corresponds to about 1800 signal events for an integrated luminosity 
of 300 ft" 1 . 

The two- lepton signal for x 2 production can also originate from the neutralino whereas 
the contamination from heavier neutralinos is very small. The total cross-section for leptonic 
X3 decays is 40 fb. Contrary to the x 2 , the two leptons from x\ originate from a real Z-boson 
and have an invariant mass equal to M z . 

For the scenario A, see Tab. 2, the production of neutralinos x 2 an d X3 nas been simulated 
with CompHEP 4.4 [45], using CTEQ6M parton distribution functions. The production 
cross-section is substantially modified by QCD corrections [48]. However, for the deter- 
mination of superpartner masses, only the kinematic properties of the decay products are 
important, which are modified relatively little by radiative corrections. For the purpose of 
this work, radiative corrections have thus been neglected. Information about superpartner 
masses can be extracted from kinematic edges in invariant mass spectra of the final state 
particles [46,49]. The distribution of the di-lepton invariant mass m« in x 2 decay has a sharp 
upper edge 

m«/,max,2 = rn^o - m^o, (26) 

see Fig. 1. The peak at mu = M z comes from the contribution of while events at lower 
invariant masses originate mainly from the x\- Assuming 300 fb -1 luminosity, a simple fit 
to the upper edge of that region gives 

m llm2 = 73.5 ± 0.5 ± 0.08 GeV, (27) 

where the first error is statistical, while the second error accounts for the systematic error 
from energy scale uncertainty in the detector (see [47] for discussion) . The error is comparable 
to what was found in [47] for the MSSM scenario (/?). 
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Figure 1: Fit to m« distribution for light neutralino production at the LHC. Backgrounds 
from Standard Model sources are not included, as they are expected to he small. 



For further studies, the decay chains involving the \z can be separated from the x° by 
applying the cut \mu — M%\ < 10 GeV on the di-lepton invariant mass. Including the jet 
from the squark decay b — * b^l gives additional information. For the decay chain with the 
X° , the invariant rriju^ distribution has an upper endpoint with 



m 



jll,max,3 



2ml 



\m 2 - m 2 - - m% - m% 



X3 



+ mtomf + mioM* + m(M, 



-(m|o - ml) ^A(m| s ,m|o,M|)] . 



(28) 



with A(a, b, c) 



b 2 + c 2 — lab — lac — 2bc. Since the mass difference between m, and 



m~ b2 is small, no experimental distinction between the two states can and needs to be made. 
Flavor-tagging of the b-jet from the sbottom decay does not improve the analysis, since the 
main background is ti. 



In a typical supersymmetry event, there are multiple jets. The jet from b — > b\^ is 
expected to be relatively hard -©rj ^ 200 GeV, but there are additional hard jets from 

the decay of the other sbottom and from gluinos, g — > bb. This introduces an irreducible 
combinatorial background. However, including that background, the characteristic edge in 
the rriju^ distribution at rriju^^ is still visible, see Fig. 2. The combinatorial background 
can be reduced by special techniques [46,47], but here we simply choose to fit it. The fit 
result is 

m jllm3 = 463.6+^ ± 2.3 GeV, (29) 
where as before the second error includes lepton and jet energy scale uncertainties. A second 
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(a) (b) 




rriju [GcV] m jU [GcV] 

Figure 2: Fits to the rriju distribution for (a) X3 an d (b) x\ production at the LHC. 
Backgrounds from Standard Model sources are not included, as they are expected to be 
small. 



edge in the rrijiig distribution is found at 

2 ^ T 2 2 4 22i22i2?\ / r2i2?\ / r2 

m i«,mi„,3 = «-2- K? m xg ~ m n ~ m x? m ?> + m n m b + m xi M z + m b M z 

x°s ' (30) 

+K° -m?)^A(m|o,m|o,M|)], 

which can be fitted in the same way as the upper end point, yielding 

m^,min,3 = 256.2+6;° ± 1.3 GeV. (31) 

In addition to studying the decay chain with the x% by requiring the invariant mass of the 
lepton pair to be sufficiently below the Z pole, vrtu < M% — 10 GeV, the decay chain with the 
X2 can be selected. Similarly to the X3 case, the mju^^ distribution has a characteristic 
endpoint at 

m 2 jllm2 = -2-04° - ™!?)K 2 - m li>- ( 32 ) 



mt 



As the X2 decays through an off-shell Z* , the mju jmax ^ distribution has no characteristic 
endpoint towards the lower end. To first approximation, the spectrum of X2 decays via an 
off-shell Z* can be thought of as superposition of Breit-Wigner line-shapes, which are close 
to Gaussian. Consequently, the upper end of the mju tmaXt 2 distribution can be approximated 
by an error function. A fit gives the rather poor result 

m jllm2 = 447l£ ± 2.3 GeV, (33) 

which is limited by statistics and the shape of the distribution near the endpoint, which is 
less steep than for the di-lepton distribution. 
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Figure 3: Correlation between m^o and m^o from LHC measurements. 



X2 



Light charginos xf can be detected in the squark decay chains by looking for a same-sign 
lepton signal originating from the processes 



pp 
pp 
pp 
pp 



99 
99 
99 
99 



bbb*b* 
bbbb - 



ttt*t* 



tttt 



bbtt W + W + XiXi 



-> bbttxtxt 

■ bbtt XiXi -> bbttW'W 
■* tt bb XiXi -> tt w _ 

itbbxtxt - ubbw + w + xixi - 



-+bbiti + i + wxiXi, 
blttrr mxixi, 
^ttUrr mxixi, 

ttbbl+l+VMXiXi, 



(34) 



see Ref. [50]. However, since besides the neutralino as the lightest supersymmetric particle, 
the neutrino in the chargino decay also escapes detection, the remaining lepton-jet invariant 
mass distributions do not allow a meaningful mass extraction. 

The measurement of the heavy neutralinos xl an d xl at the LHC is very difficult. As 
pointed out above, the appearance of a lepton pair in the neutralino decay is the best 
possibility for detection. However, due to small branching ratios of the heavy neutralinos 
into leptons, the statistics for this channel are very low. 

From the combination of the results in eqs. (27), (29), (31), and (33) one can extract the 
following absolute values for the superpartner masses, 

m-o = 33+^ 5 GeV, m^o = 106.5±^J GeV, m^o = 181±?g GeV, 



m ~ b = 499^? GeV. 



(35) 

The large errors are due to large correlations between the mass parameters, as illustrated for 
one example in Fig. 3. This can be explained by the fact that all measurements of kinematic 
endpoints in the decay distributions are closely related to mass differences, whereas no 
independent direct measurement of one of the masses, e.g. the lightest neutralino mass, is 
available. 

The analysis in this section has been performed for the specific parameter point A (see 
Tab. 2). However, most of the results are expected to be rather typical for nMSSM scenarios 
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that explain baryogenesis and dark matter. These two constraints predict that all neutralino 
and chargino states should be relatively light. In addition, our LHC analysis relies heavily on 
the gluino and a sbottom state with large left-chiral component to have masses of the order 
of a few hundred GeV. At present, the mass of the gluino is not constrained by any data, 
but if GUT-scale gaugino unification is realized, the gluino should have a mass of about 
1 TeV or less. As discussed before, the presence of light squarks are helpful in avoiding 
the suppression of chiral charges leading to the generation of the baryon asymmetry in this 
model. Following Ref. [16] we have therefore assumed that the third generation squarks have 
masses of the order of 500 GeV. While, as emphasized above, this assumption is crucial in 
the analysis of the LHC phenomenology, it plays only a minor role in the ILC related one. 
At an e + e~ collider such as the ILC, a precise measurement of the lightest neutralino mass is 
possible, so that by combing results from the LHC and ILC, the precision of the neutralino 
mass determination can be greatly improved. This was discussed already for the MSSM in 
Ref. [46] and will be analyzed for the nMSSM in the next sections. 

3.2 The nMSSM at Giga-Z 

In nMSSM scenarios that account for baryogenesis and cold dark matter, the lightest neu- 
tralino mass is typically smaller than half the Z-boson mass, m^o < M%/2. In particular 
this is the case in the reference point A (cf. Tab. 3). As a consequence, the lightest neu- 
tralino contributes to the invisible Z width. Current limits from LEP [51] already pose 
strong constraints on the nMSSM parameter space [25]. Nevertheless, these constraints can 
be improved at a future high-luminosity linear collider (ILC) running on the Z pole [52]. 
This "Giga-Z" option for the ILC could collect 50 fb _1 of integrated luminosity at y/s ~ M z , 
which is an improvement of several orders of magnitude compared to LEP. 

Since in a typical scenario like reference point A, the neutralino dark matter annihila- 
tion proceeds almost exclusively via s-channel Z-exchange, the precise determination of the 
invisible Z width opens up the opportunity to directly and model-independently determine 
the relevant ZxiXi coupling 

l^o| 2 = ^(|7V 13 | 2 -|7V 14 r) 2 . (36) 

However, the achievable accuracy at Giga-Z for r^ v is limited by systematics, and can 
reach at best 0.1% [53], which is only a factor 4 better than current LEP constraints. This 
corresponds to an error of about 0.2 x 10~ 3 for BR[Z — > XiXi]- m the case of scenario A, 
BR[Z — > XiXi] ~ 0-3 x 10~ 3 , so that the dark matter relic density can be determined from 
the invisible Z width with an uncertainty of at best 60%. Here it is assumed that the mass 
m^o is determined from some other observable with a much smaller error. 

Despite the large error, a Giga-Z analysis can confirm the existence of a neutral, (quasi-)stable, 
weakly interacting particle and a rough estimate of its coupling to the Z boson in a model- 
independent way. This would also be an interesting cross-check for more detailed model- 
dependent measurements at higher energies. 
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3.3 The nMSSM at ILC 



3.3.1 Higgs physics 

The CP-even Higgs bosons can be produced in e + e~ collisions through radiation off a Z 
boson, e + e~ — > S{Z . This process is very similar to the case of the Standard Model Higgs 
boson, except for the different Higgs- Z couplings strength as in eq. (23). The Si and S 2 
states have relatively large couplings to the Z boson, see eq. (24), and are therefore produced 
in sizeable rates. 

As mentioned above, both Si and S 2 decay predominantly invisibly into the lightest 
neutralino in scenario A. Nevertheless, the kinematic mass peaks of the Higgs bosons can 
be reconstructed from the recoil of the Z, which is cleanly characterized by the leptonic Z 
decays Z — > e + e~ and Z — > [54,55]. Due to photon radiation and detector resolution 

effects, the reconstructed mass peaks are smeared out somewhat, but their width is less than 
10 GeV, so that the two states Si and S 2 can be clearly distinguished. 

Taking into account the reduced GzzSi couplings with respect to the SM, the Higgs masses 
can be determined from the Z recoil spectrum with the precision [54] 

5M S i « 130 MeV, 5M S2 « 185 MeV. (37) 

Based on similar studies for the SM-like and invisible Higgses [55,56], the branching ratios 
of Si and S 2 can be extracted with the following accuracy, 

BR[Si -> 66] = (8 ± 0.7)%, BR[Si -> inv.] = (91 ± 3)%, 

BR[S 2 -> 66] = (2 ± 0.3)%, BR[S 2 -> inv.] = (79 ± 5)%, ' °" ! 

BR[S 2 -> W + W~] = (17 ± 1.5)%. 

Here only the statistical error is given, taking into account selection cuts to reduce the 
backgrounds. The large invisible branching ratio of both scalar states points towards a 
sizeable Higgs self-coupling A. 

The light pseudo-scalar Pi is an almost pure singlet, and the production cross-section 
e + e~ — > SiPi suppressed by more than two orders of magnitude with respect to the corre- 
sponding process in the MSSM. Furthermore, the only available final state bb+lfi is swamped 
by background from S; L Z production, so that it appears impossible to discover the Pi at the 
ILC. 

The two heavy neutral Higgs bosons S3 and P 2 as well as the charged Higgs boson are 
too heavy be produced at the ILC with center-of-mass energies up to 1 TeV. In particular 
the charged Higgs boson should always have a clearly visible decay channel into top- and 
bottom-quarks, H + — > tb. Charged Higgs production proceeds mainly through an s-channel 
photon (with only a sub-dominant contribution from s-channel Z exchange) and thus cannot 
be suppressed by modified couplings. As a consequence, the non-observation of the charged 
Higgs would clearly indicate that Ma is large, Ma > y/s/2. In this case, from the observation 
of the two light invisibly decaying Higgs states, it can be deduced that these two Higgs bosons 
must be mixtures of the SM-like Higgs and a new singlet. In other words, from the analysis 
of the Higgs bosons at the 500 GeV ILC, one can already identify an extended Higgs sector 
beyond the two doublets in the MSSM. 
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Xj = Xi 

xl 


594 32.2 



Table 7: Tree-level production cross-sections in fb at ^fs = 500 GeV with unpolarized 
beams for the reference point A (Tab. 2). 

3.3.2 Supersymmetric particles 

In scenario A, many of the neutralino and chargino states are produced with sizeable cross- 
sections at the ILC with y/s = 500 GeV, as shown in Tab. 7. The most promising production 
processes are e + e~ — > X2X.I1 X3X4) Xi~xT? xtxli with cross-sections of more than 10 fb each. 
While the lightest neutralino is not produced directly with large rates, it can nevertheless 
be studied in great detail in the decays of the heavier neutralinos and charginos. 

In the following subsections, the most promising decay channels will be analyzed in 
detail, including Standard Model and supersymmetric backgrounds. It is found that the 
two charginos and all neutralinos except the xl could be discovered and their masses and 
cross-sections measured. The discovery of a fifth neutralino state would be a smoking gun for 
a non-minimal supersymmetric model, but turns out to be very challenging experimentally. 
Nevertheless, the discovery of two neutralino states which are much lighter than the lightest 
chargino is already clear evidence for physics beyond the MSSM, where only one neutralino, 
a dominant bino-state, can be significantly lighter than the charginos. 

We simulated signal and background with the Monte-Carlo methods from [57], including 
full tree-level matrix elements and Breit-Wigner parameterizations for resonant intermediate 
particles. Initial-state radiation and beamstrahlung are always included. The processes are 
generated at the parton level, but jet energy fluctuations through parton shower and detector 
effects are parameterized by smearing functions with lepton and jet energy uncertainties 
taken from [58]. Jets overlapping within a cone with AR = v/(0x - 2 ) 2 + (77! - m) 2 < 0.3 
are combined into one jet, where fa and rji are the azimuthal angle and rapidity of jet %. 
Similarly, a lepton lying within a jet is combined into the jet. Leptons and jets outside the 
central region of the detector have a higher likelihood of mistag and get inflicted by the 
large two-photon background. Therefore leptons within an angle of | cos#| < 0.95 around 
the beam line and jets with | cos#| < 0.90 are discarded. After these simple procedures, the 
remaining isolated jets and leptons define the signature of the simulated event. 

For most processes, the signal cross-sections can be enhanced and background can be 
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P(e+)/P(e-) 


X2X4 X3X4 xtxi xtxi 


w + w~ zz tt w + w~z 


80% left / 50% right 
80% right / 50% left 


25.8 52.2 1557 51.3 
19.6 39.7 107 39.0 


24500 1020 1130 95 
770 440 500 4.5 



Table 8: Polarized tree-level production cross- sections in fb for neutralino, chargino and 
some of the largest SM background processes at y/s = 500 GeV for the reference point A 
(Tab. 2). 



reduced by a suitable choice of beam polarization. Here we assume that both the electron 
and positron beam are polarized, with polarization degrees of 80% and 50%, respectively. 
It is further assumed that 500 fb^ 1 of luminosity is spent for P(e + )/P(e~) = left/right and 
right/left each. The center-of-mass energy is always \fs = 500 GeV. The signal and main 
SM background cross-sections for polarized beams are summarized in Tab. 8. 

For all neutralino and chargino processes under study here, the main SM background 
come from double and triple gauge boson production, tt production and two-photon pro- 
cesses. Quite generally, they can be reduced by observing that the supersymmetric signal 
processes lead to large missing energy, and most of the final state particles go in the central 
detector region. In particular, by imposing a minimum value for the total transverse mo- 
mentum, p t > 12 GeV, the two-photon background is practically completely removed^ [7]. 
Furthermore, cuts on the missing energy ]fi , the polar angle of the missing momentum, 
cos# PmiBB , and of the visible momentum, cos# ptot = pi ong /Ptot are effective to reduce the 
Standard Model backgrounds. 

3.3.3 Chargino xt 

The lightest charginos almost exclusively decay via xt W ± x\i, where the W can decay 
leptonically or hadronically. Since charginos are produced in pairs, two W bosons appear 
in their decays. If both IV's decay leptonically, the rates are relatively low and one has 
to deal with a difficult background from W + W~ . The purely hadronic final state, on the 
other hand, is plagued by jet pair combination ambiguities. Therefore the best mode is: 
e + e~ — > X1X1 — > W + W~ X\ Xi ~^ 33 ^ + $ where j stands for a jet. The most important 
SM backgrounds are W + W~ , ZZ, tt and two-photon production, while supersymmetric 
backgrounds from neutralino pairs, ZSi and Higgs pairs are also taken into account. 

The SM backgrounds can be reduced by the general cuts explained in the previous section, 
p t > 12 GeV, Jj) > 100 GeV and | cos# Pmiss | < 0.8. Since the two jets in the signal originate 

§Note that the rejection of the two-photon and e ± --f background depends crucially on an excellent coverage 
of the detector at low polar angles, so that energetic fermions with low transverse momentum can be vetoed. 
The results of Rcf. [7] are based on the detector design of the TESLA study [58] , with low beam crossing 
angle, muon detectors extending to 65 mrad, and endcap calorimeters extending to 27.5 mrad. Although for 
the current ILC detector R&D several changes in the details of this setup are discussed, the planned ILC 
detector designs are expected to reach a similar photon-induced background rejection [59]. However, we also 
want to point out that the simulation of the photon-induced background in Ref. [7] with PYTHIA [60] has 
unquantified and possibly large theoretical uncertainties. 
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Figure 4: Energy distribution for the jet pair originating from the W boson in the decay 
xf — > W ± Xi after selection cuts, with a simple fit to the kinematic endpoints. 



from a W boson, the invariant jet mass is required to fulfill \rrijj — Mw| < 10 GeV, which 
removes neutralino background and is also effective on tt. The WW and tt backgrounds are 
further reduced by placing a cut on the reconstructed invariant mass of lepton and missing 
momentum, m; jmiss > 150 GeV. 

Both the signal and the main background are increased for the beam polarization combi- 
nation P(e + )/P(e _ ) = right/left, but since the signal cross-section is large, and backgrounds 
after cuts are relatively low, this polarization combination helps to increase the measurement 
precision. 

With the selection cuts listed above, about S = 105000 signal events are retained while 
only B = 30000 background events survive (mainly from Standard Model sources). The 
statistical error for total cross-section measurement is 5af x L = 0.35%. For the opposite 
polarization combination, P(e + )/P(e~) = left/right, we obtain Saf 1R = 1.3%. 

The distribution of the VF-boson energy from chargino decay, reconstructed from the 
momenta of the two jets, can be used for a chargino mass measurement. The energy spectrum 
of the W boson is almost evenly distributed with characteristic endpoints at 

' ^ m lf ~ m \\ + M w)v^T \j A(m|± , ml, , M&) (s - Ami ± ) ] . (39) 



E, 



mm. max . 9 

Ami ± 

Xi 

The energy distribution edges can be fitted by using a step function which is convoluted 
with the initial-state radiation and beamstrahlung spectrum and with the jet smearing func- 
tion. The fit function is fitted to the distribution obtained from the Monte-Carlo simulation 
with a binned x 2 fit- The fit results are (see Fig. 4) 

E min = 83.731°;°?? GeV P max = 214.8 ± 0.8 GeV. (40) 

resulting in 

m-± = 165.0 ± 0.3 GeV m^o = 33.3^ GeV. (41) 
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Figure 5: Threshold scan for chargino pair production at ILC. 

Note that only the statistical error has been given here, while the analysis of systematic 
errors at this level of precision would require a more elaborate investigation. Experience 
from the W mass measurement at LEP however shows that the systematic errors can be 
controlled to better than this level of accuracy. 

3.3.4 Threshold scan for chargino xt 

In the determination of the chargino and lightest neutralino masses from the decay distri- 
bution, the precision is limited due to substantial correlation between m-± and m^o in the 
analysis. This can be improved by using an independent measurement of m-± via a threshold 
scan. 

As an example, the measurement of the chargino pair production cross-section at six 
different center-of-mass energies below and above the nominal threshold E t ^ T = 2m-± = 330 
GeV are considered. Here it is assumed that the chargino mass is already roughly known 
from distribution measurements studied above. The measurements below the threshold allow 
to determine the background and extrapolate to values y/s > 2m -±, where the chargino ex- 
citation curve sets in. In combination with measurements above the threshold, the threshold 
energy E t h r can be precisely determined. 

It is assumed that 10 fb _1 luminosity is spent per point, amounting to total of 60 fb _1 . As 
before, the beams are polarized with P(e + )/P(e~) = right/left, and the same selection cuts 
as in the previous subsection are applied. Since the chargino mass can already be determined 
from decay distributions, this information together with eq. (39) can be used to reduce the 
background further. The result of a simulation performed with this procedure is shown in 
Fig. 5. A fit with a simple quadratic function gives a very small statistical error, 

E thT = 329.97 ± 0.1 ± 0.03 GeV, (42) 

where the second error comes from the uncertainty in the incoming beam energy, which 
is estimated to be of the order 10~ 4 [61]. This corresponds to m~± = 164.98 ± 0.05 GeV. 
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Together with the information from kinematic edges in eq. (40), the lightest neutralino mass 
is constrained to 

m -o = 33.3l££ GeV. (43) 

3.3.5 Neutralino 

Both Xz an d Xa have sizeable branching ratios into Z bosons, BR[x3 — > Zx°] = 74% and 
BR[%4 — > -Z'x?] — H%- This decay has the advantage of leading to very clear final state 
signatures, and small background contamination from SM processes involving W^-boson. 
Because of small branching ratio of the Z into leptons, only the hadronic final states are 
promising for the analysis of the neutralinos, at the cost of dealing with higher backgrounds 
and jet pair combination ambiguities. Then the signature is: e + e _ — > x^Xa ~^ ZZ XiXi ~^ 
Aj + E. 

The largest SM backgrounds arise from ZZ, W + W~ , tt and two-photon processes, but 
super symmetric backgrounds from production of other neutralino pairs, e + e~ — > xlxp with 
7^ (3,4), and chargino pairs, e + e~ — > XiXii are & l so important. The SM backgrounds 
can be suppressed by choosing P(e + )/P(e~) = left/right, while the neutralino signal remains 
sizeable for this polarization combination. With the general cuts mentioned above, p t > 12 
GeV, Ifi > 100 GeV and |pi n g /ptot| < 0.9, the backgrounds are further reduced. In addition, 
two pairs of jets have to form the invariant Z mass, |t% u - 2 — M% | < 10 GeV and \rrij 3 j A — M% \ < 
10 GeV, which removes X2X4 background and is also effective on tt. The tt background is 
reduced even further by using an anti-bottom-tag with efficiency 95% and a mistag rate of 
3% for light flavors and 25% for charm jets [62]. Finally, the chargino background is cut by 
removing events where two jets combine to give the invariant W mass, \m,j i j j — M w | < 5 
GeV. 

After these cuts, S = 400 signal events and only B = 38 background events remain, 
leading to a statistical error for total cross-section measurement of 8a^ 4 = 5.2%. 

The energy spectrum of the Z-bosons from neutralino decay has characteristic upper and 
lower endpoints both for the x° and x\i given by the expressions 

~ m l?( m l§ ~ m \l + s ) T y / A(m|o , m| ; , M|) A (m|o , rr^o , s) ) , 



(44) 



X4 



Since we consider the same decay mode for the Z bosons stemming from both neutralinos, 
a distinction between the decay products of the xi an d Xa is not possible. Therefore, the 
measured spectrum contains all four kinematic edges at the same time, which can be fitted 
with a convoluted step function as above. From the fit, see Fig. 6, one obtains 
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Figure 6: Energy distribution for the jet pair originating from Z bosons in the decays of 
the neutralinos for the process e + e~ — > x^Xa- The plots show the expected event rates after 
selection cuts, with a simple fit to the kinematic edges. 

m^o = 181.5 ± 7.6 GeV m^o — 278.0 ± 11.5 GeV, (45) 

while no good constraint on m^o is obtained. The masses of both X3 an d Xa can be deter- 
mined, but due to the low statistics, the error is relatively large. 

3.3.6 Neutralino 

The analysis of X2X4 production works similarly to the XsXa production described in the 
previous subsection. The main difference lies in the fact that the x-L due to small mass 
difference to the lightest neutralino x\i decays only through a virtual, not on-shell, Z boson. 
Again, the best statistical significance of the signal is achieved by focusing on the hadronic 
decay modes of the Z bosons, e + e~ — > X2X4 ~^ ZZ*X\Xi ~ * 4j + 1ft, and choosing the 
polarization combination P(e + )/ P(e~) = left/right. 

The backgrounds from ZZ, W + W~, tt and two-photon processes as well as X3X4 neu- 
tralino and e + e~ — > XiXi chargino production can be further reduced by the usual selection 
cuts, p t > 50 GeV, ffi > 100 GeV, |pi on g/ptot| < 0.9, \rrij i j j — M\y| < 5 GeV, and an anti-b- 
tag. Because of the large mass difference between the two produced neutralinos x\ an d X4, 
the signal typically has a relatively large total transverse momentum, so that the cut on p t 
is increased to 50 GeV here. Moreover, one pair of jets, stemming from the x°a decay, has 
to have an invariant mass equal to the Z-boson mass, while the other jet pair, associated 
with the xS], must have an invariant mass smaller then the Z-boson mass. Therefore the 
cuts \mj 1 j 2 — Mz\ < 10 GeV and Mz — iTij 3 j 4 > 10 GeV are very effective. 

These cuts reduce the background to B — 61 events, while S = 430 signal events are 
retained. Therefore the total production cross-section can be extracted with a statistical 
error of Sa^ = 5.4%. 

As before, more information can be extracted from the decay distributions. The energy 
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Figure 7: (a) Fit to the energy spectrum of 
the jet pair originating from xl decay, with 
Mi . (b ) Fit to the energy spectrum of 



the jet pair originating from xi decay, with 
<C Mz- (c) With to invariant mass spec- 
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spectra of the jet pairs stemming from the X2 an d xl decay have distinct endpoints given by 
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A-4 



(46) 



m|o(m|„-m|o + S ) T 



^Km\ v m| ; , Ml) A(m|o, m|o, s)) • 



In addition, the invariant mass of the jet pair from x\ bas a maximum value corresponding 
to the mass difference of the neutralinos x\ an d X?, 



(47) 



All these kinematic edges can be fitted with simple functions, as shown in Fig. 7, with the 
result 

m -o = 106.61$ GeV, m^o = 278.0l?| GeV, m^o = 33.31$ GeV. (48) 
As for the case of X3X4 production, the errors on the neutralino masses are relatively large. 
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3.3.7 Chargino xf 

Heavy charginos xf can be produced in association with light charginos xf j xfx^f- 
While the xf mainly decays into W bosons and the lightest neutralino xf the xf has a 
large branching ratio of 62% into ZW ± Xi v ^ a the decay chains xf ~^ ^xf ~^ ZW ± x1 an d 
xf — > W^xH — > ZW ± Xi- These decay channels for the charginos lead to three gauge bosons 
and missing energy from the neutralinos. The gauge bosons themselves can decay through 
various channels, but a good compromise between large statistics and clean leptonic final 
states is obtained by requiring one W boson to decay leptonically, and the other W and the 
Z going into a hadronic final state. The final state signature is then characterized by four 
jets, one lepton and missing energy, e + e~ — > xfxf ~^ ZW + W~XiXi ~^ 4j Z ± + If,. The 
signal can be enhanced by polarizing the beams with P(e + )/P(e") = right/left, see Tab. 8. 

The most relevant SM backgrounds are triple gauge boson processes, e + e~ — > W + W~Z 
and ti events. Production of heavy neutralino pairs can also lead to three gauge bosons in 
the final state, and thus is another important background. 

After the standard cuts p t > 12 GeV and Jf, > 100 GeV and a bottom veto, a number 
of additional selection cuts has to be applied to reduce the backgrounds. Since the SM 
background, without neutralinos in the final state, tends to give more energy to the gauge 
bosons, it can be reduced by requiring the total hadronic energy to be £had < 300 GeV. 
The invisible neutralinos in the signal also lead to a large value for the invariant mass 
of the leptonic momentum and reconstructed missing 4- momentum, m^, so that the cut 
mj^ > 150 GeV reduces the SM further. For the SM background, this invariant mass can 
also be reconstructed from the missing 3-momentum f only, by assuming that the missing 
energy originates from a neutrino, and it should be close the W mass. Therefore the cut 
— Mw| > 10 GeV is also effective on the background. Furthermore, the signal is 
characterized by a large acoplanarity between the lepton and combined jet system, so that 
the cut cos0 acOj ij > —0.7 is useful. Finally, two of the jets have to combine to the invariant 
mass of the Z boson, while the other two jets have to combine to W mass, \irij 1 j 2 — M%\ < 10 
GeV and \m,j 3 j 4 — M^\ < 10 GeV. This removes most of X2X4 background and is also effective 
on ti. 

After application of these cuts, the SM background is removed to a negligible level, while 
still a sizeable contamination of background from X3XA is left- In total B = 245 background 
events remain, compared to S — 186 events for the signal. Since the cross-section for the 
neutralino process can be measured independently, as described above, it can be subtracted, 
but the additional error from this procedure needs to be taken into account. The resulting 
expected precision for the xfxf cross-section is 5af 2 = 13%. 

For the chargino signal, the spectrum of the 4-jet invariant mass has an upper limit of 
^invj,max = m y£ ~ ui^o, which can be used to extract information about the heavy chargino 
mass. The neutralino background typically leads to slightly smaller 4-jet invariant masses, 
so that this upper edge is not contaminated. From a fit to the data, one obtains 

™ invJ , max = 287.21^ GeV, (49) 

which together with the m^o mass measurement from the analysis of XiXi production di- 
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rectly translates into 



m-± = 319.5l!;jj GeV. 



(50) 



3.3.8 Combination of sparticle measurements at ILC 

Feeding in the precise measurement of the neutralino mass from the analysis of XiXi produc- 
tion, the masses of the heavier neutralinos from X2X4 an d X3X4 production can be determined 
much more accurately, 

m^o = 106.61^ GeV, m^o = 181.5 ± 4.9 GeV, m^> = 278.0±|;| GeV. (51) 



For the lightest neutralino and the charginos, the expected errors given in the previous 
sections are not improved by combining with the other neutralino observables, so that one 
obtains 
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33.31°;* GeV, 



164.98 ± 0.05 GeV, 



X.4 



319.51^ GeV. 



(52) 



From a \ 2 fit to all mass and cross-section observables, constraints on the underlying neu- 
tralino and chargino parameters can be extracted. For completeness, we also allow a cubic 
singlet coupling k as in the NMSSM. In the nMSSM, k must be zero, but it is interesting not 
to impose this requirement a priori, but see how well it can be checked from an experimental 
analysis. The parameter k enters in the (5,5)-entry of the neutralino mass matrix, 
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(53) 



The possible measurements at the ILC analyzed here comprise mass measurements for four 
neutralino and two chargino states, as well as four cross-section measurements. They can 
be used to derive bounds on the seven unknown parameters in the neutralino and chargino 
mass matrices. Furthermore, the cross-section measurements also allow to place limits on 
the masses of the sneutrino and selectron, which appear in the t-channel of the chargino 
and neutralino production diagrams. Based on the analysis of the expected experimental 
error in the previous subsections, the following constraints on the underlying parameters are 
obtained: 



Mi = (122.5 ± 1.3) GeV, 
M 2 = (245.0 ± 0.7) GeV, 
|A| = 0.619 ±0.007, 
v s = (-384 ± 4.8) GeV, 



|k| < 2.0 GeV, 
tan/5 = 1.7 ±0.09, 
|0m| < 0.32, 



m/> o > 5 TeV, 



m e - R > 1 TeV. (54) 
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The extraction of the parameters A and v s is strongly correlated, which can be understood 
by the fact that these parameters enter in the chargino and neutralino mass matrices mainly 
through the combination \i = —Xv s . As a consequence, the effective parameter \i itself is 
determined more precisely than v s , n — (238 ±1.2) GeV. 

The results of the fit show that the sizable value of the trilinear Higgs coupling A can 
be established, which is a necessary requirement to avoid the Higgs mass bounds and allow 
a successful baryogenesis in singlet extensions of the MSSM. Furthermore, a strong upper 
bound on the value of k is obtained, which allows a distinction between the two typical types 
of singlet extensions, the NMSSM and the nMSSM. 



4 Cosmological implications 

The cosmological energy density of the main components of matter, baryons and dark matter, 
is measured with a remarkable precision [39]. In units of the critical density^ 

n B h 2 = 0.02233t°;S, 
VcDuh 2 = 0.1106lffi?, 

at 95% CL. According to the observations, the baryon density is dominated by baryons 
while anti-baryons are only secondary products from high energy processes. The source 
of this baryon-anti-baryon asymmetry and the nature of dark matter are major puzzles of 
particle and astrophysics. 

Assuming that inflation washes out any initial baryon asymmetry after the Big Bang, 
a dynamic mechanism should generate the asymmetry after inflation. Most microscopic 
mechanisms for baryogenesis fulfill the three Sakharov requirements: 

• baryon number (B) violation, 

• CP violation, and 

• departure from equilibrium. 

These conditions are satisfied in the supersymmetric extensions of the Standard Model during 
the electroweak phase transition. This is the basis for electroweak baryogenesis (EWBG). 
Baryon number violation occurs in the nMSSM due to quantum transitions between in- 
equivalent SU(2) vacua that violate the sum of baryon and lepton number B + L. When 
the electroweak phase transition is first order, bubbles of broken phase nucleate within the 
symmetric phase as the Universe cools below the critical temperature. These provide the 
necessary departure from equilibrium. To generate the observed baryon asymmetry the elec- 
troweak phase transition has to be strongly first order. For light Higgs bosons, a first order 
phase transition can be induced by the trilinear coupling between the singlet and Higgs 



^p c = 3H^/(8ttGn) where H = h x 100 km/s/Mpc is the present value of the Hubble constant, h 
0.709^032, and Gat is Newton's constant. 
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fields [16]. Such a term produces a sufficient contribution to the cubic term in the one-loop 
effective potential that is responsible for making the phase transition first order. 

As discussed in Section 2 the lightest superpartner in the nMSSM is typically the lightest 
neutralino, with a mass below about 60 GeV (for A below the perturbative bound) and a 
sizeable singlino component. Cosmologically this signals danger since a light, stable particle 
with very weak gauge couplings may have a relic abundance substantially above the observed 
amount. But if the lightest neutralino is able to annihilate sufficiently well, this particle is 
a good dark matter candidate [63]. 

For values of tan/? and A consistent with the perturbative limit, the lightest neutralino 
also acquires some higgsino component. Since its light mass allows for efficient annihila- 
tion nearby the Z° resonance, its relic density spans a wide range as the function of its 
mass. While the s-channel Z° exchange is the dominant annihilation mode, there are also 
contributions from s-channel Higgs boson exchanges generated by the trilinear singlet- Higgs- 
Higgs term of the superpotential. These contributions tend to be significant only near the 
corresponding mass poles [16]. 

When the neutralino relic density is consistent with the astrophysical constraints, and 
the model parameters with electroweak baryogenesis and the perturbative bound, the next 
lightest supersymmetric particle (NLSP) is always at least 15 percent (and frequently more 
than 25 percent) heavier than the LSP, assuring that the co-annihilation contribution is 
strongly Boltzmann-suppressed. Moreover, since the NLSP is typically a neutralino with 
large bino component, which has small couplings to the gauge bosons, the co-annihilation 
cross section between the LSP and NLSP is suppressed even further and does not play any 
role in the following discussion. 

4.1 Dark matter relic density 

The relic abundance of neutralinos is computed numerically by solving the Boltzmann equa- 
tion for the number density of the supersymmetric particles. The complex phase of A enters 
our relic density calculation directly via the couplings and indirectly through the masses of 
the neutralinos and charginos. After diagonalizing the gaugino, sfermion and Higgs mass 
matrices, we calculate the annihilation cross sections with complex couplings. In doing this, 
we follow techniques used in Refs. [64, 65] . The co-annihilation processes are checked to 
contribute insignificantly to the final result. 

After superpartners are discovered and their properties being measured at colliders one 
has to assure the consistency of the collider and astrophysical data. A crucial part of this is 
to ensure that the lightest, stable supersymmetric particle provides a reasonable amount of 
the observed cold dark matter. 

As discussed before, the LHC will restrict some of the soft supersymmetric parameters 
within certain ranges. Using these ranges, we can calculate the possible amount of neutralino 
dark matter, Qh 2 , within the given supersymmetric model. In this section we use our results 
obtained above for scenario A. To obtain an estimate of the precision the LHC can determine 
Qh 2 in the nMSSM, we randomly sampled the nMSSM parameter space in the following 
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Figure 8: Neutralino relic density as the function of the neutralino mass. Dark (light) 
blue dots represent the 1 (2) a precision of the LHC determination of Qh 2 , while gray dots 
would be allowed by LHC data, but are excluded by current low-energy and astrophysical 
bounds. Red dots show the expected ILC precision for the examined model point. The 
present WMAP and SDSS combined 2 o limits are shown by the green shaded band. The 
right frame shows the ILC scan in more details, with contours of constant values of the 
mixing parameter (|-/V 14 | 2 — |iV 13 | 2 ) indicated by the yellow lines. 



parameter region: 

< Mi < 200 GeV, 100 GeV < M 2 < 300 GeV, < |A| < 1, - vr < <p M < vr, 
-1000 GeV <v s < -100 GeV, - 1000 GeV < k < 1000 GeV, < tan/? < 30. (55) 

Additionally, for the first generation sleptons, we use the following ranges both in our LHC 
and ILC scans: 

0.5 TeV < M eR < 10 TeV, 0.5 TeV < M Ll < 10 TeV. (56) 

The results of this scan, within the parameter region that would be allowed by LHC mea- 
surements at the 1 a level, are projected on the lightest neutralino relic density vs. mass 
plane in Figure 8. 

Blue (gray) dots represent models that are (not) allowed by either the WMAP, the 
electron EDM, or WIMP direct detection limits. The dark (light) blue dots show the 1 
(2) a precision of the LHC determination of Qh 2 . All blue points satisfy all known low 
energy, collider and astrophysical constraints. Among these, the most stringent bounds 
come from the LEP and Tevatron Higgs and sparticle mass limits, WMAP, and the direct 
WIMP detection experiments. 
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Fig. 8 (left) clearly demonstrates that from LHC measurements no meaningful constraints 
on flh 2 can be derived. In contrast, as depicted in Fig. 8 (right), the ILC 1 a precision, 
indicated by the dark red dots, is comparable to the present WMAP and SDSS combined 
2 a limits, shown by the green shaded band. Even at the 2 a level, the ILC could provide 
important information. In more quantitative terms the constraints on the neutralino relic 
density from the ILC measurements are 

0.105 < nti 2 < 0.178. (2<x) (57) 

As can be seen from Fig. 8, the computed dark matter density is strongly correlated with 
the lightest neutralino mass. Therefore a precise measurement of m^o is essential for the 
accurate prediction of ^cdm- As shown in the previous chapter, for the sample scenario A 
the neutralino mass can be constrained best from chargino observables, in particular using 
a chargino threshold scan. 

Once the neutralino mass from the ILC and their relic abundance from astrophysical mea- 
surements are known with the precision indicated by Fig. 8, we can even extract information 
on the neutralino mixing matrix. This is possible because the relic density of neutralinos 
depends only on the neutralino mass and the -Z^iX? coupling, given by eq. (36). Constant 
contours of (|A r 14 | 2 — |A^ 13 | 2 ) in Fig. 8 indicate the precision at which this coupling can be 
extracted from the combination of collider and astrophysical data in the ILC era. While in 
principle it is possible to determine the neutralino mixing matrix completely from ILC data 
alone by reconstructing the neutralino mass matrix, the combination of collider and astro- 
physical data allows to extract information about the mixings without having to assume the 
particle content of the nMSSM/NMSSM. This can provide an interesting cross-check about 
the structure of the underlying model. 

The accuracy in the interplay of the collider and astrophysical measurements can be 
considered spectacular. If a consistency of the measurements as shown in Fig. 8 (right) 
were demonstrated it would be an impressive test of this model. On the other hand, it 
is worth noting that an inconsistency of the ILC and astrophysical data would pose an 
interesting, and potentially fruitful, situation. A collider indication of higher than WMAP 
or PLANCK allowed dark matter density, for example, may signal non-standard cosmology. 
Such an inconsistency could be the first sign that the temperature of the radiation dominated 
epoch before Big-Bang nucleosynthesis is not high enough for neutralinos to reach kinetic 
and chemical equilibrium. Alternatively, the entropy of matter and radiation may not be 
conserved or late production of entropy tips its balance. If instead the relic density implied 
by collider data were too small then this could be interpreted as a hint of additional dark 
matter constituents besides the neutralino. In either case, the ILC has a potential to uncover 
far reaching cosmic connections. 

4.2 Dark matter direct detection 

Another important requirement is the consistency of collider and dark matter direct detection 
data. Direct detection experiments search for weakly interacting massive particles via their 
elastic scattering off nuclei by measuring the nuclear recoil. There are numerous existing and 
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Figure 9: Left frame: Spin independent neutralino-proton elastic scattering cross sections 
as a function of the (lightest) neutralino mass. Blue and gray dots show the results of the 
LHC scan projected to this plane. Red dots correspond to the ILC scan. The regions above 
the various curves are excluded by the indicated experiments. Right frame: Spin dependent 
neutralino-proton elastic scattering cross section compared to various direct detection exper- 
iments, and to SuperK indirect limits for neutrinos stemming from neutralino annihilation 
in the Sun. The red region indicates the prediction from ILC measurements. 



future experiments engaged in this search [66]. These experiments uniformly express their 
observations in terms of the neutralino-proton scattering cross section. Most of the current 
direct detection experiments are primarily sensitive to scalar interactions of neutralinos with 
nuclei, and typically the most stringent limits are the spin independent ones. This is because 
for heavy nuclei the spin independent cross section, being proportional to the squared mass 
of the target nucleus, is highly enhanced compared to the spin dependent one. 

While both spin dependent and independent scattering rates receive contributions from s- 
channel squark exchanges, the dominant processes are t-channel Z boson and Higgs exchange 
for spin dependent and independent interactions, respectively. Due to close to resonant Z- 
exchange the spin dependent scattering cross section is enhanced in the nMSSM. However, 
as a result of the small ZxiXi coupling, the cross section will not be within reach of the next 
generation of direct detection experiments. 

Fig. 9 shows the neutralino-proton elastic scattering cross sections in the nMSSM as the 
function of the (lightest) neutralino mass. The results of our LHC (blue and gray dots) 
and ILC (red dots) scans are projected to this plane, and the regions above the various 
curves are excluded by each indicated experiment. Due to the large singlino component in 
Xi, the spin dependent neutralino-proton elastic scattering cross section is very small. The 
spin independent cross section, on the other hand, is enhanced through the sizeable scalar 
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self-coupling A, which generates the dominant contribution for the coupling of Xi to the 
CP-even Higgs bosons. As the left frame of Fig. 9 shows, the current and next generation 
direct detection experiments can already test what role the nMSSM plays in the explanation 
of dark matter, see also Ref. [67]. In particular, SuperCDMS and Xenon should be able to 
discover or rule out benchmark point A in the near future. 

The right frame of Fig. 9 shows the results of the scan for the spin dependent neutralino- 
proton elastic scattering cross section usd- The sensitivity of many spin dependent direct 
detection experiments [68], together with the SuperK limit is also shown. Contrary to the 
spin independent case, the present and near future direct detection experiments fall short 
of observing this scenario. The strongest limit on ctsd presently comes from the SuperK 
search for an excess of high energy muon neutrinos from WIMP annihilation in the Sun. 
Since, among other astrophysical quantities, this neutrino flux depends on &sd, the SuperK 
measurement is able to constrain it from above [69] . In the future, the most stringent bounds 
on dsn may also come from indirect detection experiments. Unfortunately, the sensitivity of 
the IceCUBE neutrino telescope has only sensitivity for energies above the maximum energy 
of neutrinos expected from the nMSSM scenario, so we do not expect a strong constraint 
from it. 

One should keep in mind that in general the direct detection cross sections depend on 
the CP- violating gaugino phase M , but at present there are large theoretical uncertainties 
associated with the determination of the value of 0m preferred for baryogenesis. The choice 
0m — 0.14 in the benchmark scenario A serves as a guideline, but the correct value might 
differ by a factor of two. When changing the phase to the rather extreme case 0m ~ tt, 
the spin independent cross section increases moderately by roughly 50%. Thus our main 
conclusions, that asi can be tested in the near future, while <7sd is likely out of reach of all 
planned experiments, remain the same irrespective of the value of 0m- 

The relatively small dependence of the direct detection cross section on the phase 0m of 
the gaugino mass parameters can be undertstood by observing that the lightest neutralino 
is mainly singlino and has only a small gaugino component. As pointed out in section 2.2, 
however, a relevant phase can also appear in the singlet-Higgs sector, for example in the 
parameter a\. Such a complex phase could have a much stronger impact on the direct 
detection scattering cross sections, depending on the value of arg(aA). It would be interesting 
to study this possibility and its interrelation with the generation of baryon asymmetry. We 
reserve such a study for a future investigation. 

4.3 Baryogenesis 

As shown in the previous section, the measurements of the chargino and neutralino sector 
at the LHC and the ILC provide a test of the presence of light charginos and neutralinos, 
necessary to generate the dark matter relic density. In order to probe the mechanism of 
electroweak baryogenesis with collider results, two conditions need to be tested: the type of 
the electroweak phase transition must be strongly first order, and there must be CP violating 
processes active during this phase transition. 

In our benchmark scenario, CP violation is introduced in the baryon-number generating 
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processes through light chargino currents. For this mechanism, the charginos need to be 
light enough so that they are not decoupled at the temperature of electroweak symmetry 
breaking. Using the experimental results from LHC and ILC, see eqs. (35) and (54), the 
existence of sufficiently light chargino can easily be tested. In addition, the presence of a 
complex CP- violating parameter in the chargino sector is required for baryogenesis. However, 
even with the high precision of ILC, only an upper bound on the phase 0m of the gaugino 
mass parameters can be obtained, see eq. (54). 

To test the other condition, the strength of the first order phase transition, the Higgs 
sector of the model needs to be analyzed. The strength of the phase transition can be 
calculated from the effective Higgs potential, see e.g. [16]. It depends crucially on the super- 
symmetry breaking term m 2 and a\, which are not constrained by the analysis of charginos 
and neutralinos. However, as we will show below, information about these parameters may 
be obtained by the precise determination of the CP-even Higgs boson masses, which would 
be possible at the ILC. 

From the condition of the strongly first order phase transition, one finds the following 
conditions on the parameters of the Higgs potential [16]: 



where we have introduced the quantity D for abbreviation. 

As stressed above, constraints on these parameters can be obtained from the Higgs 
masses. To relate the masses to the underlying parameters, the Higgs mixing matrices 
need to be reconstructed. Following the discussion in section 2.2, we have assumed CP- 
conservation in the Higgs sector, to that mixing occurs only between Higgs boson with the 
same CP quantum numbers. The heavy Higgs states S 3 , P 2) H* 1 with masses of the order of 
Ma are not within reach of either the LHC or ILC (the most promising discovery channels 
for heavy Higgs bosons at the LHC are suppressed due to the small value of tan /3 in our sce- 
nario). From the fact that the charged Higgs boson would not be observed at the ILC with 
y/s = 1 TeV, one could derive the limit Ma > 500 GeV. Due this bound the heavy CP-even 
Higgs boson is essentially decoupled, so that the lightest CP-even Higgs has only sizeable 
mixing with the second-lightest CP-even Higgs. Both of these masses can be measured at 
ILC. The mass matrix is 



m 2 s = -a\Viv 2 /v s - t s /v s - AV 

G {(50 GeV) 2 , (200 GeV) 2 } for perturbative A < 0.8, 



(58) 



EE 1 (A 2 t s / 

ml VA 2 /4 sin 2 2/3 + g 2 /8 cos 2 2/3 V 




(59) 




where AM| represents radiative corrections. The largest corrections stem from top-stop 
loops, but they are relatively small except for the (l,l)-entry, 
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A rough constraint on m s can be obtained by considering that one of the two eigenvalues 
of the mass matrix is always larger than the diagonal entries of the matrix, while the other 
eigenvalue is smaller than the diagonal entries. Ergo: 

M si < m s + x2y2 < M k- ( 62 ) 

Using the uncertainties for A from the fit to the ILC measurements in the neutralino sector, 
one gets 

(41 ± 3) 2 GeV 2 < ml < (114 ± l) 2 GeV 2 , (63) 

compared to the input value of m s = 106.5 GeV. Thus this very crude estimate is already 
sufficient to establish m s in the required range. 

A more detailed determination would be possible by looking at the whole matrix eq. (60), 
including full one-loop corrections [13], not only the leading term in eq. (61). The radiative 
corrections add an additional uncertainty due to the parametric dependence on the stop 
masses mi , A t and Ma- 

These parameters would need to be constrained from experiment. While in our scenario 
the stops are too heavy to be produced at a 1 TeV linear collider, one can try to search for 
them at the LHC. The study of Ref. [70] finds that a signal from decays of gluinos into stops 
can be identified with a dedicated analysis. Ref. [70] also proposes a strategy to measure 
the stop mass, although a translation to our scenario is not straightforward. Nevertheless, 
to exemplify the improvement that such a stop analysis could bring for the understanding 
of the Higgs sector, we here simply assume that the stop masses can be measured with an 
error of 5mi = 50 GeV. For the parameter A t the situation is more difficult, since it cannot 
be measured directly. However, given that in our scenario there is only a relatively small 
difference between the two stop masses, which we assume can be measured, one can infer that 
A t < 500 GeV. As far as Ma is concerned, only a lower limit of 500 GeV could be obtained 
in our scenario, as pointed out above. With these constraints, and taking into account the 
expected errors for all relevant masses and parameters, the full one-loop analysis yields 

a x = (373l£) GeV, m s = (106 ± 18) GeV, 

t s V3 = (1561 2 !) GeV , 1^1 ~ 1-0 ± °- 65 - ( 64 ) 

Note that the parameters a\ and m s can be constrained very precisely from the measurement 
of Higgs masses at the ILC. On the other hand, the necessary condition \D\ > 1 cannot be 
proven with sufficient precision, although the result in eq. 64 is consistent with this condition. 

In summary, measurements at future colliders can allow us to establish the chargino 
and Higgs mass parameters to be in the range required for electroweak baryogenesis in the 
nMSSM, but they do not seem sensitive enough to yield definitive answers to the questions 
of the first order phase transition and of the presence of additional CP violation. 

4.3.1 Electron Electric Dipole Moment 

A necessary requirement of the electroweak baryogenesis scenario is the presence of non- 
vanishing CP-violating phases in the chargino-neutralino sector. In this work, we have 
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Figure 10: Comparison of the current bound on the electron electric dipole moment with 
parameter regions allowed by expected LHC and ILC measurements for the scenario A. The 
results are given as a function of the complex phase 0m ■ 

assumed that these phases are associated with the gaugino sector of the theory. However, 
the collider, or the dark matter constraints described in the previous sections are not suffi- 
cient to determine the exact value of the CP-violating phases necessary for the generation 
of the baryon-antibaryon asymmetry. An important question is if one could obtain informa- 
tion about these phases from the measurement of, for instance, the electron electric dipole 
moment. It is advantageous to use the electron EDM since it is precisely measured, has 
relatively low theoretical uncertainties, and for the phases relevant to the model under study 
gives the strongest constraint. Since both the baryon asymmetry and the electron EDM in- 
crease with sin(0M), the electron electric dipole moment d e provides an important constraint 
on the realization of this electroweak baryogenesis scenario. 

For non-vanishing phases in the gaugino sector, the supersymmetric contribution to d e 
may become large and severe limits on the nMSSM parameter space can be obtained. Figure 
10 demonstrates that most of the LHC scan, for which 0m deviates substantially from zero 
or it, is excluded by the present 2 sigma upper limit \d e \ < 1.9 x 10~ 27 e cm. Since neither the 
LHC nor the ILC will detect the first generation sleptons if their masses are large, we allowed 
these masses to vary in the scans in a wide range: 1 < M eR < 10 TeV, 2.5 < Ml 1 < 10 TeV." 

"Barring accidental cancellations between different contributions, smaller values are only allowed by the 
present EDM limits if the phase 4>m is very small. In principle this case could not be excluded by LHC data 
alone, since neither the slepton masses nor the CP-violating phase can be well constrained. However, as can 
be seen in figure 10, even after imposing these ad-hoc bounds on the slepton masses, the LHC data would 
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For the LHC only those models survive the \d e \ limit which either have small values of 0m, 
very large values of the slepton masses, or where the one and two loop contributions to d e 
accidentally cancel. Unfortunately, since this cancellation can happen at any value of 0m, 
the EDM limit combined with the LHC data cannot shed light on the actual value of the 
phase 0m- 

New experiments have been proposed which are expected to improve the electron EDM 
limits by orders of magnitude in the next few years [71,72]. If baryogenesis is driven by a 
single gaugino phase of the nMSSM such as studied in this work, then a non-vanishing value 
of d e will probably be measured by the time of the ILC operation as scenario A suggests. 
This can even happen if the first generation sleptons are very heavy, as shown by the case of 
the input model A, where the first generation sleptons are fixed at O(10TeV). If an electron 
EDM is measured, the indirect ILC determination of d e will also be an essential cross check 
confirming the nMSSM. According to Figure 10, the ILC data alone will be able to constrain 
|0m| below about 0.3 and \d e \ below about 5 x 10~ 28 ecm. The correlation between \d e \ and 
|0m| is strong in the relevant 0m ~ 0.1 region. Utilizing this, even a 50% measurement of 
| d e | combined with the ILC data will be able to constrain M in a ±0.15 window. 

4.3.2 Comment about parameter space for CP violation 

In the previous chapters we have made the simplifying assumption that the only CP- violating 
parameter beyond the CKM matrix is a common phase of the gaugino mass parameters 
M = 0Mi = 0m 2 — 0m 3 - However, as pointed out in chapter 2, the general nMSSM 
with minimal flavor violation can have 13 complex CP- violating parameters in total, four 
of which can be relevant for electroweak baryogenesis. These four independent parameters 
can be chosen to be Mi, Mi, a\ and A t . For illustration, here we list the different scenarios 
that are possible if these phases are allowed to be non-zero. A more detailed investigation of 
these scenarios would be beyond the scope of this paper and has to be performed elsewhere. 

Gaugino phases without gaugino unification: In this case, the baryogenesis could be 
driven by bino currents, which are sensitive to the phase of Mi, while the phase of M 2 could 
be zero. While, depending on the mass pattern, this still could allow the generation of a 
sufficiently large baryon asymmetry [36], it would reduce the contribution to the EDMs. As 
a result, lower values of the slepton masses of the order of a few TeV would be allowed. 
The collider sensitivity to 0m x would be similarly weak as for the common gaugino phase 
0m discussed above, since one would rely on the same observables in the neutralino sector. 
Contrary to the MSSM, a phase in Mi will have no sizeable effect on the dark matter 
annihilation, since in the nMSSM the LSP is expected to be mainly singlino. 

Only CP-violating phase in ay. As shown in Ref. [27], a sizeable phase in the Higgs 
sector could succesfully explain electroweak baryogenesis. With the freedom of R- and 
PQ-transformations, this phase could be attributed alternatively to a\, t s or m.12, without 
changing the physics. In this situation, the EDM constraints from one-loop diagrams will 

allow nMSSM scenarios for which the value of the electron EDM varies over several orders of magnitude. 
Therefore the parameter space scan presented here should give a representative picture of the information 
that may be obtained after running the LHC. 
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be much weaker than for a gaugino phase, since one-loop diagrams with a Higgs boson 
are suppressed by an additional small Yukawa coupling. As a result, in this scenario, light 
slepton masses of a few hundred GeV would be allowed, which could be in reach of the LHC or 
ILC. This feature could be used to experimentally discriminate this scenario from a scenario 
with non-zero gaugino phases. For instance, if the LHC or ILC would discover selectrons, 
the present electron EDM limit would rule out phases in the gaugino sector that are large 
enough to explain electroweak baryogenesis (except for the possiblity of large cancellations 
between two contributions to the EDM). A CP-violating phase in the Higgs sector could 
lead to observable effects in Higgs production processes and decays. For example, mixing 
between CP-even and CP-odd Higgs states would enable the possibility to produce a third 
light Higgs boson ^ in the process e + e~ — > <piZ. A more conclusive answer to this question 
would require further study. As pointed out in section 4.2, a complex phase of a x can also 
strongly affect the dark matter scattering cross sections and thus could relax the constraints 
coming from dark matter direct detection experiments. 

Only CP-violating phase in A t : If the only non-zero phase (in addition to the CKM 
matrix) was in the parameter A t , it is unlikely that a sufficiently large baryon asymmetry 
could be generated, as experience from the MSSM shows [34-37]. It is however possible that 
a non-zero phase of A t can exist in conjunction with other complex nMSSM parameters. 

CP-violating phases in several parameters: When considering to most general scenario 
with CP-violating phases in multiple parameters, the interpretation of the EDM limits be- 
comes more involved. From our analysis, the prospects for disentangling the complex pa- 
rameters through collider data do not seem promising. 

5 Conclusions 

In this article, we have presented the phenomenological properties of the nMSSM, assuming 
the parameters to be close to the ones that lead to a solution of the dark matter and baryo- 
genesis problems. A light neutralino and chargino spectrum appears under these conditions, 
that can be probed at the LHC and the ILC. At the LHC, however the detection of these 
weakly interacting particles becomes difficult, unless they are produced from the cascade 
decay of strongly interacting particles. In order to study the LHC phenomenology we have 
therefore assumed that the gluino is within the reach of the LHC, as suggested by gaugino 
mass unification, and that the third generation squarks are light, which is helpful in avoiding 
the suppression of the chiral charges necessary for baryogenesis. We have shown that, under 
these conditions a relatively good determination of the LSP and other neutralino masses 
may be obtained. However, the accuracy of these measurements is not sufficient to allow 
the computation of the neutralino relic density with any meaningful precision. A definitive 
probe of this model can only be done at the ILC. At the ILC, for a representative point, we 
have shown that both chargino and four of the five neutralino masses may be determined. 
The sparticle mass pattern leads to a good discrimination of this model from the MSSM, 
and even a distinction between the nMSSM and the NMSSM, which has a different singlet 
sector, is possible. Moreover, the mass measurements, together with the production cross 
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section information, lead to a good discrimination of this model from the MSSM, and al- 
low to compute reliably the annihilation cross section of the LSP, and thus to check the 
agreement needed for a successful explanation of the dark matter relic density. 

The requirement of a sufficiently strong first order electroweak phase transition also trans- 
lates in this model into the presence of a light singlet scalar, which mixes after electroweak 
symmetry breaking with the MSSM Higgs doublets. Consequently, two light CP-even Higgs 
scalars appear in the spectrum which couple with reduced couplings to the weak gauge 
bosons, and which decay predominantly into the LSP. These CP-even Higgs bosons may be 
searched for at the LHC and the ILC in the invisible decay channel. At the LHC, however, it 
will be difficult to discern between one or more invisibly decaying Higgs bosons. A definitive 
scrutiny of this question may only be performed at the ILC, at which a good determination 
of the Higgs boson masses may be established. A light CP-odd Higgs boson also appears 
in the spectrum, of difficult detectability due to its strongly reduced couplings to the SM 
fermions. Nevertheless, the pattern of the two CP-even Higgs boson masses allows to distin- 
guish the nMSSM Higgs sector from the MSSM. In addition, the mass measurements provide 
information about the parameters of the model relevant for a strong first order electroweak 
phase transition. Note that for this kind of analysis, we have shown that the large radiative 
corrections in the Higgs sector need to be under control, requiring some information about 
the superpartners of the top quark needs from measurements at LHC or ILC. 

Finally, we have also investigated the constraints coming from the current direct detec- 
tion searches for dark matter, as well as the prospects of a successful observation of the dark 
matter candidate of this model at future experiments. We have shown that the predomi- 
nantly singlet component of the LSP makes its direct detection easier than in the MSSM for 
the spin independent channel, but more difficult for the spin depedent channel. For the spin 
independent case, the current and next generation experiments should be able to definitely 
probe this scenario. Similarly, assuming that the dominant phases are in the chargino sector, 
we have also investigated the bounds coming from the electron electric dipole moment. We 
have shown that, due to the small values of tan/5 necessary to realize this model, the EDM 
constraints become weaker than in the MSSM, and, again, this model may only be probed 
by next generation experiments. If the dominant phases are in the Higgs sector instead, 
the EDM bounds are even weaker. Such complex phases may also affect the predicted dark 
matter direct detection cross sections significantly. In this case additional information about 
CP- violation might be obtainable from collider data in the Higgs sector. 

In summary, the nMSSM provides an exciting framework for addressing the problems 
of baryogenesis and dark matter. The properties of this model and role it plays in solving 
these problems could be probed with high precision at the next generation of laboratory 
experiments, and would allow us to make connections between laboratory and astrophysical 
observations at a new level of insight. 
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